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Abstract. In this paper we obtain the wave equation modeling the nematic liquid- 
crystals in three space dimensions and study the lifespan of classical solution to 
Cauchy problem. The almost global existence to classical solution for small initial 
data will be presented. 
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1. Introduction 

1.1. Physical background. Liquid crystals are a state of matter that have properties 
between those of a conventional liquid and those of a solid crystal that are optically 
anisotropic, even when they are at rest. In this paper, we shall study a variational wave 
system from liquid crystal. 

At the beginning, we introduce the physical derivation by the least action principle for 
the variational wave equations. The mean orientation of the long molecules in a nematic 
liquid crystal is described by a unit vectors n = n(t, x) £ § 2 , the unit sphere, where 
(t, x) £ R + x R 3 , which could be modeled by below Euler-Largrangian equations derived 
from the least action principle (TJ [21] , 

, , 8 fW(n,Vn) , 

(1.1) «n tt + jun t -| = An, nn=l, 

on 

where the well-known Oseen-Franck potential energy density W is given by 

(1.2) W (n, Vn) = ±a(V • n) 2 + \f3 (n • V x n) 2 + ± 7 |n x (V x n)| 2 . 

Here the positive constants a, /3 and 7 are elastic constants of the liquid crystal, corre- 
sponding to splay, twist, and bend, respectively. The Lagrangian multiplier A is deter- 
mined by the constraint n ■ n = 1. The inertia and viscous coefficients k and \x are two 
non-negative constants, respectively. 

There are many studies on the constrained elliptic system of equations for n, and on the 
parabolic flow associated with it, where we refer the reader to see [21 El EJ [T3l EI]. In 
particular, the one-constant elliptic model leads to the equation for harmonic maps taking 
values in the two-sphere [TJ [5]. One parabolic system of equations for n coupling the 
compressible Navier-Stokes equations in one space dimension (1-D) has been considered 
by 0. 

However, for the complete hyperbolic system (|l.lj) . only the 1-D case has been sys- 
tematically studied, where the main difficulty is the possible gradient blowup. For the 
1-D Cauchy problem of (|1.1[) with \i = and k = 1, which describes the model with 
viscous effects neglected, the global weak existence and singularity formation have been 
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extensively studied. In fact, for this extreme case, one example with smooth initial data 
and singularity formation (gradient blowup) in finite time has been provided in |10j , and 
the global weak existence results have been provides in (3] [4j [28j |29l [30] [31] . In a recent 
paper [5], the singularity formation (gradient blowup) and global weak existence for the 
1-D Cauchy problem for the complete system (|l.lj) (k and /i are non-negative constants) 
with a — 7 have been established. This result shows that for 1-D Cauchy problem we 
should essentially expect the similar existence and regularity results for the complete 
system (jl.lj) and the extreme case with /i = and K = 1. 

Very recently, Li, Witt & Yin [17) have proved the singularity formation and studied 
the life-span of the classical solution for a 2-D axisymmetric model. However, the well- 
posedness for (|1.1[) in 2-D and 3-D is still wide open. In this paper, we consider the 
life-span of the classical solution for the 3-D Cauchy problem of (jl.ip with small initial 
data. Note the initial data in the 1-D singularity formation examples in [5] [10] have large 
C^-norm. 

1.2. The Variational Wave Equation. In this paper, we focus on the 3-D Cauchy 
problem of (jl.ip with fx = and k = 1 and planner deformation: 

(1.3) J {{d t n) 2 -W{n,Vn)} dxdt = Q, |n| = 1, 
where the planer deformation means that 

(1.4) n = (rai, n 2 , n 3 ) — (cosu(t, x), sinit(t, x), 0) 
with the angel u = u(t, x), (t : x) e M+ x M 3 . 

Now, the variation (|1.3|) can be reformulated as 

(1.5) '{^(n f ) 2 -W(n,Vn)}dxdt = 0, 
which gives rise to the Euler-Lagrange equation 

3 

(1.6) d tt ni + d ni W(ii,Vn) -J2dxi(dd jni W(n,Vn)) =0. 

i=i 

Multiplying three equations in (|1.6p by — sinu and cosu, then adding them up, we can 
get a variational wave equation, by (|1.4I) and (|1.2I) , 

(1.7) u u - cli^d^u - c\(u) dl 2 u - f3dl 3 u + (a - 7) sin2u d Xl d X2 u = F(u,du), 
where 

Ci(u) 2 = a sin 2 u + 7 cos 2 u, c 2 (u) = a cos 2 u + 7 sin 2 u, 

and 

(1.8) F(u,du) = -(a -7)sin2u [(9iu) 2 - (d 2 u) 2 } - (a - 7) cos 2u d x u d 2 u. 

Instead of the direct calculation from (jl.6p . another way to derive (|1.7p is from the 
variational principle 

S- f\(d t n) 2 -W(n,\7n)\dxdt = 0, Inl = 1. 
ou J 

where in the derivation we need to use (|1.4I) and 

(1.9) W — ^-(cosud X2 u — sinud Xl u) 2 + ^r(d X3 u) 2 + ^-(cosud Xl u + sinud X2 u) 2 . 
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Since a, (3, and 7 are all positive constants, we can rewrite (jl.7l) as the following equation 
by scalling 

3 

(1.10) u tt - aij{u(t,x))didjU = F(u,du(t,x)), x G R 3 , t G R + , 

where 9i = ^7, (i = 1, 2, 3) and 

(c\(u) i(a — 7)sin2w 

i(a-7)sin2u <%(u) 

1 

It is obvious that there exist constants Mq,M±, such that 

< M < cf(u), c 2 2 {u) < Mi. 

Furthermore, A is a positive definite symmetric matrix. 

In our paper, we shall study the lifespan for the variational wave equation (jl.7l) in 
three dimensional space with small and smooth initial data 

(1.12) u(0,x)=e f(x), u t (0,x) =e g(x), 

where eo is a constant small enough, and (/, g) G C°°(R 3 ) with compact support. We 
will obtain the following almost global existence result: 

Theorem 1.1. The equations (ll.7[) and (|1.12[) or equivalently f| 1 . 10[) and (|1.12[) admit 
a unique solution in C°°([0, T*), R 3 ) with 

(1.13) T*>C*exp— , 

£0 

where C is a constant independent ofso- 

We have already known from [5] [10] that there are 1-D gradient blowup examples when 
the initial data have large C 1 norm, even when the amplitude oscillation of the initial 
data is small. We note that the 1-D examples are also special examples in 3-D with 
initial data restricted in one space dimension. More precisely, for the blowup examples 
in [SJUnj, in which they also consider the planar transformation n = (cosw, sinu, 0), the 
initial data uq = u(0,x) satisfy ||uo||l°= = 0(e), H^oHc 1 = 0(1) and the blow-up time is 
at O(l), where e > is an arbitrarily small number. Clearly, these singularity formation 
results of [5) [10] do not conflict with our result. On the contrary, combining these two 
pieces of information from O [10] and this paper, we have a fairly complete picture of the 
smooth classical solutions for the 3-D liquid crystal equation (|1.10[) . 

We remark for the axisymmetric case. In |17j , the authors have proved the singularity 
formation and the lifespan for the axisymmetric case in 2-D with small data. Due to 
the significant difference between 2-D and 3-D for the wave equation, for the 3-D ax- 
isymmetric case, whether the solution of the equations (|1.7p with small initial data form 
singularity or not is still open, although the nonlinear structure of our model is similar 
to the one in Lindblad [23]. We leave this case to be concerned in the future paper. 

To prove Theorem 11.11 we use the generalized energy method which is popular in 
analyzing wave equations. There are large amount of literatures focusing on this topic. 

See mm US] US [2D1 [23 [23] and [TS] f ° r m ° re details ' In Lei-Lin-Zhou QJ5], the authors 
have presented a global existence for Faddeev model in both 2-D and 3-D cases, where 
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the equations have a null condition structure. However, the equation (I1.10[) does not 
have the null condition structure, so we can not get the same results. 

This paper is organized as following: in the next section, we recall some well-known 
results for the wave equation, and in Section 3, we shall prove our main results. 



2. Preliminary 

Without loss of generality, in the following we assume the positive constants a and 7 
satisfy 

a < 7. 

The proof for the case a > 7 is very similar. Then, by a schedule of re-scaling, we can 
rewrite the p,10[) in the following wave equation 

3 



(2.1) 

here and hereafter 



Uu— Yl aij(u(t,x))didjU = F(u,du(x,t)), igR 3 , t > 0, 



t = 0, u = ef, u t = eg, 



(2.2) U =d 2 t -Y J d: 
and 



3 

,2 

Xi ' 

i=l 



(2.3) A = (a l:) )i<ij<3 = 7jCosin2w c sin 2 u 



Co cos 2 u icosin2u 
sin 2u Co sin u 




(2.4) 



is positive semi-definite with Co = > 0. 

In the beginning, we introduce the following notations. Write 

x = t, 
d = -d/d t , 
di = d/d Xi , 

d = (d Q ,d 1 ,d 2 ,d 3 ) = (-dt,d 1 ,d 2 ,d 3 ). 
Define the first order differential operators as 

(2.5) ft = (n a b)o<o,b<3, 
with 

(2.6) fl ab = x a d b - x b d a , (0 < a, 6 < 3), 
and the scaling operator as 

3 

(2.7) L = ^ / x i d i . 

i=0 

Denote the vector field 

(2.8) r = (d,n,L). 

It is well known that the above operators T have possessed a good commutation with the 
wave operator □ = <9 2 - A (See also [T41IT8] ): 
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Lemma 2.1. For multi-indices <;, we have 

(2.9) [n,r«]= ^ r ' D > 

|c|<IS|-l 

where [ , ] stands for poisson bracket and are constants. 
For any integer N > 0, denote 

(2.10) Ht,-)\\ r ,N^= £ H r Mv)IU- 

| S |<JV 

Then, we have the following decay estimate 

Lemma 2.2. Suppose that h = h(t, x) with (£, x) £ R + x M™ is a function with compact 
support in the variable x for any fixed t > 0. For any integer N > 0,we have 

(2.11) \\h(t, -)||r,iv,oo < C(l + i)-i(l + |f - |a;||)-5||ft(t, -)[[r,iv + [ ? ]+i,p, 
where p > 1, C is a constant. 

Remark 2.1. TTie proof of Lemma 2.2 can be found in [21] . 

In order to prove Theorem II. 1[ we also need the following estimate for a linear wave 
equation (See [20"1 122] ) : 

Lemma 2.3. Let h(t, x) and w(t,x) be the functions with support < t + p} in the 

variable x for any fixed t > 0, if for all norms on the right hand side of below inequality 
are bounded, then holds 

(2.12) \\hdw\\ L 2<c p \\vh\\ L 2 £ ||r Q HU- 

q|<1 



3. THE ALMOST GLOBAL EXISTENCE 

In this section, we will prove Theorem ll.il Denote 

(3.1) E S , T = {v = v(x,t)\ sup ||<9u||r>,2 < £, with v(0,x) = e f, v t (0,x) = e g}. 

0<t<T 

where s,Eo and e are given positive constants with s > 9 and eo, e small enough. We 
shall prove our main result by finding a fixed point in E St T- First, we give the linearized 
equation of (I1.10[) as the following case: 

.Ou- J2 aij{v{t,x))didju = F(v,dv(x,t)), xeM 3 , t>0, 

(3.2) { ij= i 

t = Q,u = e f,u t = e g, 

where v € -B s ,t is any given function. Therefore, the equation (|3.2I) is a linear wave 
equation, which admits a unique global solution. We define a map 

Mv = u. 

It is sufficient to prove the map M is a contract map. To do this, we apply the operator 
r fe on both side of (|3.2[) (here and hereafter k is a multi- indices with |fc| = s). By Lemma 
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2.1, we have 

3 

(3.3) ar k u- a ij (v(t,x))d l d j r k 

i,3 = l 



\a\<\k\ i,j=l 

where [ , ] is the poisson bracket. Multiplying d t T k u and taking integral on R 3 for (|3.3 
we get 

(3 - 4) ^(ll rfeau lli 2 (^)+ E f ^{<t 1 x))a i r k i M d j r k u i dx 

C ^ I r a F(v,dv(x,t))d t T a udx 

\a\<\k\ Jr3 
3 

V / [a ij (v(t,x)d i d j ,T k ]ud t T k udx 

frit Jr? 



+ \ [ d t a lj (v(t,x))d l T k u l d j T k u l dx 

R 3 



djdij(v(t, x))diT k Uid t T k Uidx, 
R 3 



2 
1 

~ 2 

where Ck a are constants. Then, we get 

3 

(3.5) W h du\?mF&) + J2 / a ij (v(t ) x))d i T k u i d j r k u i dx 

<Ce + c( V / / r Q F(v,di<x,i))<9 t r a ucte* 
V |«|<|fc| - 70 Jr3 

+ H I I [a^(v(t,x),r k }d l d J ud t T k u l dxdt 



R 3 JO 
1 ' « 



dtOij(v{t, x))dil UidjI uidxdt 



, dja i j(v(t,x))d i r k Uid t r k Uidxdt) 
2 Ji?3 Jo J 



=Ce 



o 

By the chain rule and lemma 2.2, we get 

(3.6) III + IV <C\\du\\l s2 V / 1 1 ^Oij ( v (t , a:) ) 1 1 i oo d* 

i' A ' 

ft 



<C||^||L,2 / ||5«||i«A<Ce||fl«||^ i2 / <r>- 1 dr, 
Jo Jo 



here and hereafter, we use the notation 

< x >= \f\ + x 2 . 
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Recalling s > 9 and \k\ = s, we have, 
3 ,t 



(3.7) // <C V f f V I* (a tf («(t, x))d i:j u)dT k udxdt 

^E/V E E (r^(^,x))^r^ 



+ |%r a ur^a y (u(t,a;))|)|9r fc u|da;di 
<C / ||v|| r)[ i*i ])00 ||9u|lr,|jfe|,2 + ll u llr,[^i],ooll 5u llr,|fe|,2||3f || r ,|fc|,2rfT 

•J 

<Ce(|0u||r,.,a + 1) ||3u||i>,2 / <r>-' dr, 



as well as 

(3.8) I<cf[ (||r a sin2w|| 2 ;oo + ||r a ft;|| i -) 

° M<[^] 

• Y, (I|r^(sin2^)|| i2 + \\T^dvf\\ L ,)]dt-\\du\\r, ]k]a 

\(3\<\k\ 

<cj^ (\\v\\ ri[i]+hoo + 0(\\vf vm+lt j \ ||ft;||r,|fc|,2dt ■ ||0u||r,|k|,2 
<Ce 2 \\du\\ r , s ,2. 

where we used the Taylor's formula. Recalling that is positive semi-definite, then 
by we S et 

(3.9) WMr,.* < Ceo + Ce\\du\\l.» ln(l + t), 



where C and C are constants. Then from (|3.9p . let T < Tq < e^c 7 ! and we take the 
initial data small enough such that 

Ce < |£ < 1, 

then 

Mv = uG E s ,t- 

One still needs to prove M is a contraction map. To do this, let Vi,V2 G £<s,t, and then 

Mvi = Ui 6 E St T, i = l,2. 

Subsequently, we need to estimate \\dui — 9tt2||r,s,2- For simplicity, we write U = ui — u%, 
and correspondingly, V = v\ — V2, then from equation (12. ip . we have 

(3.10) E (<Hi(«i(*.aO)fl|« 

ui - aij(v 2 (t, x))d l} u2 
i,j=i ^ ' 
-I- F(«i,<9i>i) - F(v 2 ,dv 2 ) 

= E ( °>ij{ v i{^ x ))d%jU ) + E ( _ «y( u 2)]9ijW 2 

i,i=l V - ' i,j=l V - 

+ F(t>i,3wi) - F(v 2 ,dv 2 ). 
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Then by the standard energy estimates, we get 

3 



(3.11) \\dU\\ r , s ,2+Yl I a ij (vi(t,x))d i r a Ud j r a Udx 



1,3 
3 ,t 



\a\<\k\ i,j=l 
3 r t 



< Yl E / / \{^ 3 {vi)d rjl T a ]Ud t Y k U\dxdt 



JR 3 



+ E I I \da ij (v 1 )\\dT k U\ 2 dxdt 
+ E f [ \r k ([a lJ (v 1 )~a lJ (v 2 )}d lJ u 2 )d t T k U\dxdt 

+ 11 \T k (F(vi,dvi) -F(v 2 ,dv 2 ))d t T k U\dxdt 
=1 + 11 + III + IV, 

where |fc| = s. Recalling that vi,v 2 6 E s ,t and using Lemma 12. 2[ we can do similar 
estimates as (f3T4l) - (f3T8l) 



n-t 3 

(3.12) I + II <C\\dU\\ 2 rs2 V \\da ij {v 1 {t t x))\\ LB odt 

+ C J2 fl E E (|r a M«i(^))^T^| 

lJ=1 JO JR* M < [i kt± ] \P\<\k\-l 

+ \dijT a UT^a i j(v 1 (t,x))\)\dT k U\dxdt 
<C"\\dU\\ 2 w \n(l + t). 

where we used Taylor's formula to deal with dij(vi), and hereafter, C" is a uniform 
constant. Doing the expansion 

2 sin(t>i + v 2 ) sinV 2 cos(wi + v 2 ) sin V 

(3.13) (dij(vi) - aij(v 2 )) = | 2cos(^i + v 2 ) sinV sin(«i + ua) sin V 



we can deal with 77/ as following 

(3.14) in<cj2 f I EE (\rMv&*)w p «*\ 

i,3=l J ° JR3 |a|<[JM]|0|<|*| V 



+ |%r Q u 2 r /3 sm(V(t,x))\ ) \3T k U\dxdt 



<C J ||9F|| r![ |k| ])0o ||9ti 2 ||r,|ife|,2 + IKII r)[ i*L ])00 l|9V r ||r,|fc|,2dr||5y|| r) | fc | i2 
<C"e{\\dV\\l St2 + \\dUf rtSt2 ) In(l+i), 
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where we use Lemma l2.3l to get the second inequality. Similarly, we can deal with IV as 

(3.15) IV<CJ2 I I Yl J2 (\F a (cos2v 1 {t,x)dv 1 )dr f} V\ 

+ \dT a VrP ( cos 2v 1 dv 1 ) \ + \T a (dv 2 cos 2(v x + v 2 j)T f} (dv 2 cos V) \ 
+ \T a (dv 2 cosV)T f3 (dv 2 cos(wi + v 2 ))\^\dT k U\dxdt 

<Cfcln(l +t)(l + f o (\W\\l^ +lloo + II%j>l +1]i0 > 

+ C J\ < r >~ 1 (1 + e < r >^ + < r || dV\\l lkl . 2 )dr^j ||9F||| )|ib|)2 

<C" £ ln(l + t)\\dV\\l^ + C" <t ||SF|| r , s>2 . 
Then from ([3JT]l - ([37TB]) . we have 

3 

(3.16) ||^||raa+y] / Otf(wi(*,a;))Sir a D'5 J r a crda! 

" : 7/l " 

<C"eln(l + t)\\dV\\l s , 2 + C'eHl + t)\\dU\\l s , 2 + C" <t > l \\dV\\^ 2 . 
Let T <Ti < e^u , then there holds 

3 

(3-17) \\dU\\l 8t2 + [ ^M(t^))W a Ud 3 T a Udx<\\\dV\\l^ 2 . 

Combining (I3.9[) and (|3.17l) . take T < mm{T ,Ti}, we get the desired result. 
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